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Abstract 

We study the bispectrum of the curvature perturbation on uniform 
energy density hypersurfaces in models of inflation with two scalar 
fields evolving simultaneously. In the case of a separable potential, it 
is possible to compute the curvature perturbation up to second order 
in the perturbations, generated on large scales due to the presence of 
non-adiabatic perturbations, by employing the 5./V-formalism, in the 
slow-roll approximation. In this case, we provide an analytic formula 
for the nonlinear parameter /nl- We apply this formula to double 
inflation with two massive fields, showing that it does not generate 
significant non-Gaussianity; the nonlinear parameter at the end of in- 
flation is slow-roll suppressed. Finally, we develop a numerical method 
for generic two-field models of inflation, which allows us to go beyond 
the slow-roll approximation and confirms our analytic results for dou- 
ble inflation. 
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1 Introduction 



A key prediction of a period of inflation in the very early universe, is the 
generation of a spectrum of primordial perturbations. Such perturbations 
naturally arise from the zero point vacuum fluctuations in quantum fields, 
which are stretched to arbitrarily large scales during inflation. The distribu- 
tion of the primordial density perturbations thus provides an important test 
of any inflation model. In particular, slow-roll models of inflation generically 
predict an almost scale- invariant, almost Gaussian distribution of primordial 
density perturbations [HI21- 

Typically one calculates the power spectrum of density perturbations and 
the relative amplitude of gravitational waves. However, there is limited in- 
formation in the power spectrum over the restricted range of scales where 
the primordial power spectrum can be reliably inferred from observations. 
As a result, there has been growing interest in calculating the distribution 
of primordial perturbations from inflation, considering not only the power 
spectrum but also the bispectrum and other measures of possible deviations 
from purely Gaussian distribution as a possible discriminant between differ- 
ent models jHlHj- For instance, it is known that the size of the bispectrum 
during single-field, slow-roll inflation is related to the spectral tilt of the 
power spectrum and is thus constrained to be small j^]. On the other hand, 
inflationary models with higher derivative operators, such as ghost infla- 
tion or inflation based on the Dirac-Born-Infeld action, can produce higher 
non-Gaussianity jEHSEHS], possibly detectable. Furthermore, non-adiabatic 
perturbations produced during multi-field inflation can certainly generate de- 
tectable non-Gaussianity in the density field after inflation, in models such 
as the curvaton scenario [TUJ (see also [TT] ) or reheating [T2l IT3] ITil ITS]. 

What is less clear is whether nonlinear large-scale evolution of pertur- 
bations during inflation is capable of producing significant non-Gaussianity. 
Such a study requires a consistent treatment of nonlinear perturbations, not 
just in the matter fields but also in the gravitational field. Recently a num- 
ber of authors have developed gauge-invariant descriptions of the nonlinear 
curvature perturbation on large scales [TBI 113 CHI UJH I2H1 ETJ • 

Rigopoulos et al. [22] have used numerical simulations of the stochastic 
dynamics on large scales and found significant non-Gaussianity even in simple 
two-field models. It is this question which we wish to address in this paper 
by presenting analytical and numerical estimates of the non-Gaussianity in 
two-field inflation models. We use the 5 N- formalism, as advocated by Lyth 
and Rodriguez [21]], to calculate the evolution of the curvature perturbation 
after Hubble exit. We find that in the case of a separable potential, it is 
possible to construct an analytical formula to express the bispectrum of the 
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curvature perturbation in terms of the potentials and slow-roll parameters 
of the two fields. Furthermore, we also use numerical solutions to go beyond 
the slow-roll approximation after Hubble exit, confirming our analytic results. 
In the case of double inflation with two massive fields, we find no significant 
non-Gaussianity produced, which confirms a previous discussion by Alabidi 
and Lyth in Ref. [23]. 

Note that the multi-field scenario considered here is different from that 
considered in ^TJ E3], where only one field is contributing to the energy 
density during inflation. 

In Sec. 2 we review the two- and three-point statistics of field perturba- 
tions during slow-roll inflation, and how in the <5iV-formalism these can be 
related to the two- and three-point statistics of the curvature perturbation 
on uniform density hypersurfaces in arbitrary multi-field models. In Sec. 3 
we show how one can actually calculate the evolution of the primordial power 
spectrum and bispectrum using the slow-roll approximation in simple two- 
field models with a separable potential, during inflation. In Sec. 4 we discuss 
the issue of how the end of inflation can affect the curvature perturbation and 
the non-Gaussianity. In Sec. 5 we present both an analytical and numerical 
study of non-Gaussianity in the double quadratic inflation model studied by 
Rigopoulos et al. [22], going beyond the slow-roll approximation. We present 
our conclusions in Sec. 6. 



2 Non-Gaussian perturbations in multi-field 
inflation 

Here we review some of the results concerning perturbations and non-Gaus- 
sianities in multi-field inflationary models. For simplicity, we will consider 
a model of inflation driven by a set of minimally coupled scalar fields with 
canonical kinetic terms and potential W, described by the action 



S = / d*x^g 



- v^d^i + w(ipx, ¥2, ■ ■ •) 
1 i 



(i) 



where m| = 8ttG is the reduced squared Planck mass. 

We consider scalar perturbations in a quasi-homogeneous flat Friedmann- 
Lemaitre- Robertson- Walker spacetime with scale factor a(t) and perturbed 
metric 

ds 2 = -(1 + 2A)dt 2 + 2B ti dtdx i + a(t) 2 [(1 - 2^)<J - + 2E tij ] dx { dx j . (2) 

Primordial cosmological perturbations are usually expressed in terms of 
the curvature perturbation on uniform energy density hypersurfaces, denoted 
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by C, defined in !2o] in linear perturbation theory, and generalized at 
higher order in the perturbations in [TBI El EE]- (See also [JjJJ |27J EH] for 
comparison between different variables at second order.) This quantity is 
widely used, especially because it is conserved, on large scales, for adiabatic 
perturbations [29J. Here we want to compute the three-point correlation 
function of ( and hence we need to define ( only up to second order in the 
perturbations, 

C = _^ _ _ E 5 p+ i#p + - 2 5 P 5 P + i f ? ) 5p 2 , (3) 
P P P 2 2pVP/ 

where p is the background energy density, 5p its perturbation, and H = a/a 
is the Hubble rate, while a dot denotes the derivative with respect to cosmic 
time t. 

According to the so called <5iV-forHialism [3UJ HHl E2| , C> evaluated at some 
time t = t c , is equivalent, on large scales, to the perturbation of the number 
of e-foldings N(t c , £*, x) from an initial flat hypersurface at t — i*, to a final 
comoving - or, equivalently, uniform density - hypersurface at t = t c . Hence, 
one has, on large scales, 

C(t c , x) ~ 5N(t c , u, x) = Af(t c , t„ x) - N(t c , Q, (4) 

where 

N(t e ,t m )=J Hdt (5) 

is the unperturbed value of N '. Equation (J3J) simply follows from the def- 
inition of Af as the volume expansion rate of the t = const hypersurface, 
integrated along the integral curve of the unit vector orthogonal to the t = 
const hypersurface 31J. This definition is not restricted to linear theory 
but holds also at second or higher order in the cosmological perturbations 

jaaiaaiEEEiEi!. 

One can then take t* as the time, during inflation, when the relevant 
perturbation scales exited the Hubble radius, k = aH, and t c as some time 
> £* during or after inflation. Then the number of e-foldings M can be viewed 
as a function of the field configuration <pj(i*,x) on the flat hypersurface at 
t = t* and of the time t c . If one splits the scalar fields in a background 
value and a perturbation, <p/(t,x) = ipi(t) + 8ipi(t,x), SN(t c , £*, x) can be 
expanded in series of the initial field perturbations 5<pi(t*, x). Retaining only 
terms up to second order, one obtains 

5N(t c , U, x) = N ,M + \H NuS^M, (6) 
/ i j 
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where 

are the first and second derivatives of the unperturbed number of e-foldings 
N(t c ,t*), with respect to the unperturbed values of the fields at Hubble 
crossing. Note that in general M depends on the fields, fi(t), and their first 
time derivatives, <j>i{t). However, if slow-roll conditions 

~ -Wj, at t = U (8) 

are satisfied at Hubble exit, then Af depends only on the field values [3*T] . 

In [20] it was shown that Eq. (Jlj with © can be used to compute ( up 
to second order in the perturbations in multi-field models of inflation. In 
particular, the 57V"-formalism is equivalent to integrate the evolution of ( on 
super-Hubble scales from Hubble exit until t c . 



2.1 Two-point statistics 

The power spectrum of the curvature perturbation (, V^, is defined as 

2tt 2 

(C kl Ck 2 ) = (27r) 3 5( 3 )(k 1 + k 2 )—V c (h). (9) 

Let us now take the scalar field perturbations as uncorrelated stochas- 
tic variables at early times, with scale invariant spectrum of massless scalar 
fields in de-Sitter space. Their two-point correlation functions hence satisfy 

WO = (2vr) 3 5 /J 5( 3 )(k 1 + k^P^), V.ik) = £|, (10) 

where if* is evaluated at Hubble exit, k = aH. From Eqs. and (fTUj) . 
making use of Eqs. (jlj) and (jBJ), one obtains 

n = E^*- (ii) 

Another interesting observable that can be derived from the two-point 
statistics is the scalar spectral index of V^, defined as 

n ^ l = TTnT- (12) 

For multi-field inflation, the expression for the scalar spectral index has been 
given, for instance, in p?T) IT]. At lowest order in slow- roll one has [2] 



d\nV c d\nV c 1 d\nV c 

2. — — 

dink ~ dN H dt 



(13) 
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The first equality follows from taking k at Hubble crossing, k = aH oc 
HexpN, while the second is a consequence of the definition of N, Eq. @. 
Making use of the expression for the scalar power spectrum given in fTTf . 
and of 

one thus obtains 

where e is the well known slow-roll parameter defined as 

H , s 

^~Jp- (16) 
This expression can be shown to be equivalent to that given in [HI]. Indeed, 
the particular combination of second derivatives, ipjNjj, that appears in the 
spectral index (fT5j) can be eliminated in favor of second derivatives of the 
potential using the slow-roll approximation to give 

" t " 1 ^' 26 "i4E^ + VT K N% ■ (17) 
2.2 Three-point statistics 

Let us now discuss the three-point statistics of the curvature perturbations. 
The bispectrum of the curvature perturbation ( is defined as 

(C kl Ck 2 Ck 3 > = (2n)W 3 \J2 k i) B dkn fe, h)- (18) 

% 

Observational limits are usually put on the so called nonlinear parameter 
/nl- We use here the definition of fc-dependent /nl given in [S], 1 

6 f - n ^ B < (10) 
~ 5 /nl = E^fl^Tf " (19) 

One can compute the bispectrum using Eqs. (jlj and For the three- 
point correlation function of £, this yields 

(c kl c k2 Ck 3 ) = E NjN^N^lMMi) + 

UK 

- £ NjNjN KL (5p{ 1 8vl 2 (6<p K *6<p L ) k3 ) + perms, 

1 UKL 

(20) 

1 Note that this choice of sign for /nl differs from that used by The minimal 

detectable /nl with an ideal CMB experiment is slightly larger than unity \WM I5T] . 
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where a star denotes the convolution and we have neglected correlation func- 
tions higher than the four-point (see [35J|3I]). 

The first line represents the contribution from the three-point correlation 
functions of the fields. For purely Gaussian fields, this vanishes. However, 
Seery and Lidsey [HE] have found that during slow-roll inflation the field 
three-point correlation functions do not vanish and read 

(s<pLMMs) = M^E^^* 2 E f^-M (21) 

u ' • perms ^ n " L V 



with 



k 2 k 2 1 k 2 k 2 



M = M(h, k 2 , hs) = -hk 2 2 - 4-p + -kf + -j^(k 2 - h) , (22) 

where kt = ki + k 2 + k 3 , and it is assumed that all the ki are of the same order 
of magnitude so that they cross the Hubble radius approximately at the same 
time. The function Ai, which has been written in a slightly different form 
from that in Ref. [S^, parameterizes the /c-dependence of the three-point 
functions. The sum is over all simultaneous rearrangements of the indices 
J, J, and K, and the momenta hi, k 2 , and k 3 in JH, such that the relative 
position of the ki is respected (3*Hj . 

The sum over the permutations of M. over all the ki reads 

J 7 (ki,k 2 ,k 3 ) = E M{h,k 2 ,k 3 ) 

perms 

-ill^ + i^M-lZki], (23) 



and one can sum over the permutations and evaluate the first line of Eq. (}2*U[) 
using Eq. pTf and 



E NjNjN k E <PihKM{hik 2 ,h) = -Hj2N}F(ki,k2,h), (24) 

UK perms 1 

where we have used J2i Njtpi = —H. 

To evaluate the second line of Eq. ([20)1 . as in [33] we assume that the 
connected part of the four-point correlation functions is negligible and we 
make use of Wick's theorem to reduce the four-point functions to products 
of two-point functions. Hence, we can finally rewrite the bispectrum using 
its definition (fT%|) and Eq. (|2T)|) . This yields, after few manipulations, 

„,, , , , ,, 4-r.2E.fc? ( - l ? T.i,,NjNjNjj \ 
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To derive this expression we have used Eq. (|TT|) to replace Vl with V%. 

One can relate the bispectrum to the nonlinear parameter /nl by using 
Eq. (HI, 

' 5 /NL " 2^ (1 + /} + (E^iVp ' (26) 

where 

/ = /(*i.k.*s) = -1-2^5 ( 27 ) 

is a function of the shape of the momentum triangle with the range of values 
< / < | [E] , and we have used Eq. (jlljl in the first term on the right hand 
side of Eq. f!26|) . The lower bound on / is obtained in the geometrical limit 
when one of the three k^s is much smaller than the other two, e.g., k\ <C 
k 2 ~ k 3 , in which case /nl becomes independent of k jSlEH], while the upper 
bound is obtained in the equilateral triangle configuration, k\ « A; 2 ~ &3. 

The first term on the right hand side of Eq. (|2fi|) is momentum dependent 
and comes from the three-point correlation functions of the fields, Eq. (|21j) . 
computed by quantizing the perturbations inside the Hubble radius during 
inflation j^SlEZI- Notice that V* can be related to the amplitude of gravita- 
tional waves. Introducing the ratio between tensor to scalar modes, 

8V 

r = (28) 

this term can be written as 

~^ = T6 (1 + /) ' (29) 

and is constrained by observations to be small, r/16 <C 1. This expression 
simplifies the bound on this term found in Ref. [37]. 
The second term |2H] . 

_ 6 (4) = EijNjNjNjj 

5 /NL " {Y.kN%Y ' m 

is momentum independent and local in real space, because it is due to the 
evolution of nonlinearities outside the Hubble radius during inflation. 

The nonlinear parameter /nl can only be larger than unity if the mo- 
mentum independent term is large. We will devote the next section to 
compute this quantity in two-field models. 
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3 Two-field inflation with separable potential 

We will consider now two scalar fields (pi = <f> and ipi = Xi described by the 
action 



S = - rj fj ^V=9 [^0^0 + \d» X d»X + W{4>, X ) 



(31) 



We assume that the potential of the fields is separable into the sum of two 
potentials each of which is dependent on only one of the two fields, 

W{cf>,x) = U(<P) + V( X ). (32) 

3.1 Slow-roll dynamics of background fields 

The Klein-Gordon equations for the background fields read 

4> + 3H<p + U' = 0, (33) 

X + 3H X + V' = } (34) 

where 

*'% « 35 > 

The unperturbed Friedmann equations read 

H = 



^(0 2 + X 2 )- 07) 



2m P 

Inflation takes place when e = —H/H 2 < 1. Here it is assumed that the 
slow-roll conditions are satisfied for both fields during inflation. In this case, 
the Klein- Gordon and the first Friedmann equations reduce to 

3H<j) ~ -U', 3H X ^ -V, 3rripH 2 ~ W. (38) 

It is hence convenient to extend the definition of slow-roll parameter e for 
one field and define f3E3 

with 

e = + e x . (40) 
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One can also define the two slow-roll parameters 





(41) 



Note that throughout this paper, for simplicity, we will assume U' > and 
V > so that we may eliminate first derivatives of the potential in favor of 
the slow-roll parameters Vr and \feX. 

We will use the slow-roll equations (|3*8|) to write the number of e-foldings 
(0) during inflation as [3UJ [T] 



There is a crucial difference between single field inflation and inflation 
with many fields [SB]- In single field inflation the slow-roll solution forms a 
one-dimensional phase space. This means that once the inflationary attrac- 
tor has been reached, there is a unique trajectory. In particular, the end 
of inflation takes place at a fixed value of the inflaton field which in turn 
corresponds to a fixed energy density. However, if two fields are present, 
the phase-space becomes two-dimensional and there is an infinite number of 
possible classical trajectories in field space. The values of the two fields at 
the end of inflation will in general depend on the choice of trajectory. In- 
terestingly, for the potential ([3^1. under the slow- roll conditions (|38|). there 
exists a dimensionless integral of motion C, using which one can label each 
slow- roll classical trajectory, 



The number of e-foldings iV in Eq. (j4*2*|) characterizes the evolution along a 
given trajectory, while the quantity C allows us to parameterize motion off 
the classical trajectory, and it will turn out to be very useful in computing the 
curvature perturbation during inflation. In this respect, the study presented 
here is very similar to the one discussed in Ref . [3S1 , in the case of a potential 
that is a separable product rather than a separable sum. 

3.2 Perturbed expansion during slow-roll 

In order to calculate £, given in Eq. (p£J), we need to calculate the perturbed 
expansion due to the field quantum fluctuations on an initial spatially flat 
slice up to a final comoving or uniform density hypersurface at t c . To calculate 
the power spectrum of £ (fTTj) we require only the first derivative of the 




(42) 




(43) 
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expansion with respect to the initial field values. But to calculate the scalar 
spectral index (JT5j) or the three-point correlation function (J2U|) we need to 
compute the perturbation in the number of e-foldings expanded up to second 
order in the initial field perturbations. We will thus proceed as follows. We 
will first compute the first derivative of N(t c ,t*) with respect to the fields, 
i.e., N t $ and N tX , by differentiating N(t c , t*) in Eq. (}4*2j) in d<fi* and dx*- With 
these two first derivatives we will be able to compute the second derivatives 
Njj. 

We first note that each of the integrals in Eq. (|4*2|) depends upon both 
0* and x*- F° r instance, the value of the integral over depends upon x* 
through the dependence of the limit C upon C (</>*, x*), t ne integral of motion 
labelling the classical trajectory. In other words one has 



dN = —r 

mi, 



mi 



U 

¥' 

v_ 

V' 



dxc 



QXc 
dx* 



v 

V' 



d X * VU' 



u_ 

u 



dx*. 



(44) 



In order to compute the derivatives of the number of e-foldings with 
respect the initial fields 0* and x*, one needs to compute dxc/°X*i dXc/d4>*, 
etc. Since the value of C and Xc on a given classical trajectory is a function 
of the conserved quantity C, one has 



dXc 



dXc 
dC 



dC , 


dC 


ha d 

0(p* 


4>* + rv 

dx* 


' dC 






OX* 



dx* 



dx* 



(45) 



Making use of Eq. 



one can easily compute 
OC mr, dC 



nip 



mi 



dx* V* 



(46) 



We now require that t = t c coincides with a p = const hypersurface, which 
during slow-roll is given by 



Ufa) + V{ X c 
Differentiating this condition yields 



const . 



dC c dC c 



0. 



(47) 



(4f 
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Furthermore, differentiating Eq. 
using Eq. (}4*%j) . one finally obtains 



with respect to the trajectory C and 



2 "Vc 

mp dc 

2 ^Xc 

Thus, substituting Eqs. 



- + - 



v' (— + — 

c I K' 2 £/' 2 



(49) 



and (jm in Eq. (gSJ), one finds 



W c I e 



dXc 



W*e c 



W*e c 



1/2 



dXc 



WceJ i < 
W*e c 



<A 1/2 



X\ 1/2 
c 



(50) 



where e is the overall slow-roll parameter given in Eq. ()16j) and we have 
written the derivatives of the potentials that appear in Eqs. (}4*6|) and (J4*9|) in 
terms of the slow- roll parameters defined in Eq. (J3H|) . 

One can now evaluate the derivatives of the number of e-foldings with 
respect to the initial fields 0* and x*i using in Eq. the expressions (joHj) . 
For the first derivative this yields 



dN 



m P - 



1 U* + Z c 



nip 



dN 

dx* 



'2e 
1 V t 



where 



(V4-U c eX)/e c . 



(51) 
(52) 

(53) 



To compute the second derivatives we differentiate Eqs. (|5Tj) and (|52~|) 
with respect to 0* and x*- This yields 



ni 



mp 



d 2 N 
d 2 N 



dxl 



rrir 



d 2 N 



V* K - Z c 
2e* W* 
nip dZ c 



nip 



dZ r 



W^2ef U( P* 
nip dZ c 

W*V2eJ d X * 
nip dZ c 



Wj2e\ 



dx* 



W*V2e 



x 



(54) 

(55) 
(56) 
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In single field inflation, perturbations on large scales are adiabatic and 
hence ( remains constant during inflation . In two-field models, however, 
the large-scale perturbations are not in general adiabatic and ( evolves dur- 
ing inflation. In the previous expressions, Eqs. (JHTj) . (f52~J) and (jniHSHl), the 
function Z c takes account of this evolution by expressing how the local expan- 
sion depends upon the local field values on the uniform density hypersurface 
during inflation. 

3.3 Curvature perturbation and non-Gaussianity dur- 
ing inflation 

We now compute the curvature perturbation and the non-Gaussianity on a 
comoving or uniform energy density hypersurface during slow-roll inflation, 
as given in Eq. (j3J). From the first derivatives of N, Eqs. ()51|) and (p)2"j) . one 
can derive an expression for the power spectrum of the curvature perturbation 
. It can be simplified by introducing two dimensionless variables, related 
to the values of the potentials and first slow-roll parameters at and t c , 

U* + Z c V* — Z c 

u ^ w, ' BS -wr- (57) 

The power spectrum then reads 

(u 2 v 2 \ , . 

Note that this can never be less than the power spectrum derived by con- 
sidering only adiabatic perturbations restricted to the background trajectory 
(at fixed C). In multi- field inflation there is an additional contribution to 
the power spectrum due to non- adiabatic perturbations at Hubble exit. Thus 
one has V( > V^\ c , where Vc\ c can be obtained by taking the limit t c — > £* in 
Eq. (|EBJ), since for adiabatic perturbations the curvature perturbation does 
not change after Hubble crossing. In this limit 

ef 6* 
u — > — , v — > — , for t c — > t*, (59) 

and one finds 

To compute the scalar spectral index with Eq. (15) and the nonlinear 
parameter /nl with Eq. (26) one needs the second derivatives iV/j, given in 
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Eqs. ([54Hod]) . In the expression for the scalar spectral index the derivatives 
of Z c with respect to the fields contained in these equations mutually cancel. 
Thus, by using Eq. (flHjl . one finds 



u 



n f - 1 = -2e* - 4^ ^f— , . (61) 

4 ^ 

Note that in the limit t c — > t*, using Eq. (|55|). one finds the spectral index 
for purely adiabatic perturbations, 

n c|o - 1 = -6e, + 2rT, (62) 

where we have defined jSHl 

^ CTCT _ ( £< y + e x^x)/ e- ( 6 3) 

This represents the effective mass of adiabatic fluctuations tangent to the 
background trajectory. 

The momentum dependent nonlinear parameter, f^, defined in Eq. (j29jl . 
takes a very simple form, derived from Eq. (|6()|) . 

-f/$ = e.^(l + /)<<*(! + /), (64) 

where / is the momentum dependence, defined in Eq. (12 7|) . 

In order to give an expression for one needs dZ c /d(J)* and dZ c /dx*, 
which can be computed by employing Eq. (|50j). One finds 

Ve*7r— = -V 6 *^ — = — ^*"4> (65) 
00* c>x* m P 



where we have defined 



e 

with 



^-Zr,— . (66) 



^ = (^ + e ^x)/ e _ ( 67 ) 

This represents the effective mass of isocurvature fluctuations orthogonal to 
the background trajectory. 
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We are thus able to give an analytic expression for f^l defined in Eq. (fHU|) . 
One finds 

v> ■ jM-Hfr-^o+te-*)^ (68) 

- tJnl - 1 7 T! • l b °J 

This is the exact expression for the amplitude of the nonlinear parameter 
for the curvature perturbation ( during slow-roll inflation in an arbitrary 
two-field inflation model with separable potential, and represents one of the 
main results of this paper. 

As for the power spectrum and spectral index, we can take the limit of 
t c — > £*, which gives the nonlinear parameter for purely adiabatic perturba- 
tions in the two-field case. Using Eq. (JHHJ) in Eqs. (fMJ) and (jEHJ), yields 

-j?/NL| c = e*(l + /) + 2e*-rr (69) 

which coincides with the single field case result EE]. Note that, using 
Eqs. (JH2|) and (|2Kj) for the adiabatic case, this can be expressed as a consis- 
tency relation between observables [5], 

-|/NL| = ~(ne| -l) + ^/, (70) 

and is thus constrained to be small. 

Although the general expression (p3%|) for two fields, which allows for non- 
adiabatic perturbations, is rather involved, we can qualitatively discuss the 
order of magnitude we expect for f^. Since max(u, v) < 1 and u + v = 1, 
the denominator in Eq. (jHSj) is of order e~ 2 , where we use e to denote generic 
first-order slow-roll parameters, e or rj. The first two terms in the numerator 
of Eq. are of order e~ l , so their contribution to is of order e*. Only 
the third term, which is of order Ae~ 2 , leads to a contribution that is not 
automatically slow-roll suppressed. 

Although t] ss may become larger than unity during inflation, the prefac- 
tor in front of the parenthesis in Eq. (jfifij) can be very small or vanishing, 
suppressing the contribution of A to the nonlinear parameter. If this is not 
the case, since A does not appear in the expression for the spectral index 
rif, Eq. (jHIl), for r\f je c 3> 1 one may have models with large / NL and quasi- 
scale-invariant spectral index, corresponding to large deviations from the 
consistency relation for purely adiabatic perturbations, Eq. fTOjl . We leave 
to future work the investigation of models where A, and thus /nl, are large. 
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4 Perturbations after inflation 



So far we have considered only the curvature perturbation during inflation. 
In order to relate our calculations to observables we need to calculate ( after 
inflation when the universe is radiation dominated, on a uniform energy 
density hypersurface for t = t c > t e , where t e denotes the end of inflation. In 
this case t c defines a uniform radiation density hypersurface. 

In single field inflation large-scale perturbations are adiabatic and thus the 
curvature perturbation on uniform density hypersurfaces remains constant 
both during and after inflation, independent of the detailed physics occurring 
at the end of inflation. In two-field models we need to take account of how 
the local expansion depends upon the local field values both during inflation 
and at the end of inflation in order to calculate ( some time after inflation. 

If one of the two fields - for instance <fi - has stabilized before the end of 
inflation, so that the end of inflation is dominated by a single field \i then Z c 
becomes constant, and the power spectrum, the scalar spectral index and the 
non-linear parameters after the end of inflation are simply given by Eqs. (J58j) . 
fH}, tH) and (JHEJ) with Z C = U C = const. Note that in this case A = at 
the end of inflation and the nonlinear parameter f^ will be small. 

If both fields are "active" at the end of inflation, and inflation takes place 
on a hypersurface q((p e , Xe) = const at t = t e , undefined - i.e., not necessarily 
a uniform energy density, one can separate <5iV(t c , t*) into the sum of two 
pieces, the first from t — i* to the end of inflation t — t e , the second from 
t — t e to a uniform density hypersurface after inflation t = t c , 

6N(t e , £*) = 5N(t e , Q + 6N(t c , t e ). (71) 

The first piece, 5N(t e , t*), can be computed by a similar calculation to that 
of Sec. 13.21 It is possible to show that 5N expanded at second order is found 
by replacing Z c in Eqs. (|3tT)I . (1521 and (jHlHHfij) by 




Then one needs to compute 5N(t c ,t e ). 

In scenarios such as the curvaton or modulated reheating scenarios, it 
is assumed that the end of inflation hypersurface is effectively unperturbed, 
5N(t e ,t*) ~ 0, and that the primordial density perturbation originates from 
isocurvature field perturbations, 5s, during inflation, which introduce a per- 
turbation 5N(t c , t e ) oc 5s only after or at the end of inflation. 

An alternative possibility is to assume that inflation ends due to a sudden 
instability triggered by some function of the fields reaching a critical value 
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|40[ ET] • This is what happens in hybrid inflation |l3*] El] where the 
false vacuum state is destabilized when the inflaton field, 0, reaches a critical 
value, = (f) e . If we assume instantaneous reheating of the universe one has 

1331 

5N(t e ,t e ,x) = N '(P(*c,x)) - N '(p(*e,x)) 



—Sp+-—5p 2 
dp 2 <ip 2 




pW 
p/ p 2 



(73) 



We assume that the energy density is conserved at the end of inflation, so 
that the background value of the energy density in this expression is p(t e ) = 
U((j>e) + V{Xe)- Furthermore, if the equation of state becomes radiation-like, 
so that p = —4Hp, we then have 



SN{t c ,t e 



Sp 
P 



2 p 2 



(74) 



Note that the numerical coefficient that appears in front of the brackets in 
this expression is dependent upon the equation of state after inflation has 
ended and would be 1/3(1 + w) for an equation of state P = wp. 

In a hybrid-type limit, where the false vacuum dominates the self- interac- 
tion energy of the slowly-rolling fields, we can take the effective potential 
to be almost completely flat (consistent with the slow-roll approximation) 
so that [5p/p] e is negligible and 5N(t c ,t e ) 5N(t e ,t*). The primordial 
curvature perturbation ( is then given by the curvature of the end of inflation 
hypersurface, ((t c ) = 6N(t e ,t#). 



5 Double quadratic inflation 

To be more specific and give a quantitative estimate of the non-Gaussianity, 
it is instructive to consider the case of massive fields with potential given by 
Eq. with gg HH] 

U = \m\^\ V = \m\ X 2 ■ (T5) 

These scalar fields thus have no explicit interaction, but can interact gravi- 
tationally during inflation. 
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5.1 Slow-roll analysis 

We can apply our earlier analysis to calculate the curvature perturbation 
during slow-roll inflation in this model. 
The integrals (021) and (g3J) yield 

m 2 P N(t c , Q = + xl) ~\(<t>l + xt) , (76) 

mfC(t c , Q = m-/ In - m" 2 In , (77) 

where, as for N, we have fixed the limits of integration in the definition 
of C in Eq. (|43J) to run from to t c . In the standard treatment, one can 
parameterize the slow-roll trajectories of the scalar fields in polar coordinates 

m 

X = 2m P v / iVsin6>, = 2m P VN cos 6. (78) 



The advantage of doing so is that the angular variable 9 can be related to 
the number of e-foldings by the expression |1H| EB1 ITT] 

N = N n 1 j (79) 

where we have defined the ratio between the masses of the fields as 

R = m x /m ( f ) . (80) 

Thus, Eq. (J79|) implies that an inflationary model is completely specified by 
giving R and the values of the fields <fi and x ~ or equivalently iV and 9 
- at some given time. Once specified an initial condition, Eq. (|79J) can be 
used to evolve the background fields or equivalently the slow-roll parameters. 
Then one can use Eqs. (JSEJ), flU), (IMj) and (fRHjl. to evaluate completely 
analytically the power spectrum V^, the scalar spectral index and the 
nonlinear parameter /nl for double inflation. 

Before studying in detail these equations for a specific case, one can es- 
timate the amount of non-Gaussianity generally produced during and after 
inflation. During inflation, as explained in Sec. \'<\.'<\\ /nl can be large only 
if, in Eq. (68), the contribution proportional to the quantity A defined in 
Eq. (}6T)j) is large. In double inflation this contribution can become temporar- 
ily large, at most of order unity, at the turn of the trajectory in field space. 
However, after inflation, the fields have settled to their minima, which im- 
plies Z c = and also ^4 = 0. In this case the dependence on the field values 
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at t — t c in Eq. (joTj) disappears and we can simply rewrite Eqs. (|55jl. (pTTj) . 
and (J58j) by using the potential (J75)l and Eq. (|75jl with C = Xc = 0, as 




N(t c ,U) 



(81) 



47r 2 mp 



1 



(82) 



— 1 



— 2e* — 



N(t e ,UY 
— [2 + f(k 1 ,k 2 ,h 




2N(t c ,U) 



1 



)], 



(83) 



where N(t c ,t*) is the number of e-foldings between Hubble exit and the 
moment when 4> c = Xc = 0, which roughly coincides with the end of inflation. 
These results coincide with the estimate of Ref. [2H]. We conclude that double 
inflation cannot produce large non-Gaussianity - i.e., |/nl| "C 1- these being 
suppressed by the slow-roll conditions on the fields at Hubble exit. We now 
turn to a more detailed analysis of Eqs. (JH%|) . (jEH), (jHljl and (JHBJ), and compare 
with a numerical results for the evolution on large scales without assuming 
slow-roll. 

5.2 Numerical analysis 

The 5iV-formalism assumes that the fields are in slow- roll at Hubble exit, 
Eq. (jHJ) (see [IB] for a development of the <5iV-formalism applicable to more 
general situations). Furthermore, in deriving Eqs. (|58|). (jfiTj) . (JMJ), (J55|) and 
the background evolution relation ([79)1 . we have assumed slow-roll all along 
the inflationary evolution, from i* to t c . 

While still assuming that slow-roll conditions are satisfied at Hubble exit, 
one can go beyond the slow-roll approximation by numerically solving the full 
background evolution equations for the fields, Eqs. (|3^j) and (jSljl . together 
with the Friedmann equation (|36|). and computing the evolution of the scale 
factor a. One can then evaluate the expansion iV = In a as a function of 
the field initial values (0*, x*), up to some final time t c which may be some 
time during or after slow-roll inflation. Then, one can numerically calculate 
the first and second partial derivatives of N(t c ,t*) with respect to 0* and 
X* by the finite- difference method jlHl, i-e., by computing N for different 
values close to (0*, x*)- Indeed, this provides an efficient method to compute 
the non-Gaussianity numerically for any two-field model, also when the field 
potential is not separable. 

We study, as an example, the specific case considered by Rigopoulos et 
al. in [22]: mass ratio R — 1/9 and initial condition <p = x = 13mp. We 
assume that inflation ends on a hypersurface e e = 1 and that Hubble exit 
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Figure 1: Examples of trajectories in field space are shown from Hubble exit, 
starting 60 e-foldings before the end of inflation at e e = 1, for R = 1/9. The 
thick (orange) line represents the trajectory starting from = x — 13mp, 
and shown from Hubble exit, (0* = 8.2, x* — 12.9), to the end of inflation, 
(0 e = 0.0, Xe = 1-4). The grey shading represents the space of all possible 
trajectories. 
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Figure 2: The values of the fields (solid line) and x (dashed red line) during 
the inflationary trajectory of Fig. ^ are shown as a function of the Hubble 
rate H, from iV ~ 42 to iV ~ 4 e-foldings from the end of inflation. Note 
that time increases from right to left. 




Figure 3: The energy densities of the fields (solid line) and p x (dashed 
red line) normalized to (mpm^) 2 , are shown as for Fig. |21 
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Figure 4: The power spectrum V( of the large-scale uniform density pertur- 
bation C during the inflationary trajectory of Fig. ^ is shown as a function 
of the Hubble rate H, from N ~ 42 to iV ~ 4 e-foldings from the end of 
inflation (mp = 1). The solid and the dashed lines represent the analytic 
and numerical calculations, respectively. 
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Figure 5: The spectral index is shown as for Fig. 01 
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Figure 6: The nonlinear parameter / NL , in the limit where f(k\, k 2 , £3) = 
(k\ <C k 2 ~ fe) is shown as for Fig. EJ The (turquoise) dotted line represents 

— (6/5) /nl /(! + /) while the (blue) dashed-dotted line represents — (6/5) /^lj 
both computed analytically. 



takes place at i = i*, corresponding to 60 e-foldings before the end of infla- 
tion. One finds that (0 e = 0.0, Xe = 1-4), corresponding to 9 e = n/2, and 
that (0* = 8.2, x* = 12.9), corresponding to 6>* = 0.327T. 

The trajectory in field space is shown in Fig.^ One can see that assuming 
slow-roll from £* to t e is well justified. Indeed, the turn of trajectory in field 
space, corresponding to the moment when <fr exits slow-roll, takes place near 
= 0, when the energy density of </> has already become negligible. This 
can also be checked in Figs. |2] and El where the values of the fields and their 
energy densities are shown, respectively, at the turn of the trajectory in field 
space, from N ~ 42 to iV ~ 4 e-foldings from the end of inflation. 

Now we compute and /nl = /nl+/nl analytically, using Eqs. (|55j) . 

(jHIJ), (EU) and (JHEl), and numerically, using the procedure explained above. 
In Figs. EHElwe compare the analytic and numerical calculation of the power 
spectrum the scalar spectral index n^, the nonlinear parameter /nl, as 
a function of the Hubble rate during inflation. The two calculations agree 
remarkably well, within an accuracy that depends on the numerical precision. 

As expected, when the heavy field <ft dominates the total energy density, 
for H/m^ > 0.8, the uniform density curvature perturbation ( remains con- 
stant. It starts growing when H drops to 0.55 < H/m^ < 0.8, corresponding 
to both the light and heavy fields contributing to the total energy density. 
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During this phase, the large scale total entropy perturbation is non- vanishing, 
and sources the evolution of ( 38 . After few oscillations, the heavy field en- 
ergy density is red-shifted and the light field dominates the universe while ( 
becomes constant in time again. Consequently, the power spectrum and the 
spectral index grow monotonically, following the growth of 

The evolution of the nonlinear parameter — /nl> however, is non- mono- 
tonic. It grows sharply during the intermediate phase 0.55 < H/m^ < 0.8, 
corresponding to the heavy field leaving slow-roll, but it then decreases. This 
growth corresponds to the growth of A defined in Eq. (IfifiJ) . 

One can compute the value of /nl at the end of inflation, for e e ~ 1, in the 
limit where f(k u k 2 , k 3 ) = (h < k 2 » h). One finds -(6/5)/j$ = 0.008 
and — (6/5)/nl = 0.008 which coincide with the analytical result of Eq. ([8*3^1. 
— (6/5) /nl = 0.016 ~ 1/60. We conclude that /nl in this model is much 
smaller that unity. 

Note that, as shown in Fig. 01 the curvature perturbation becomes con- 
stant after the energy density of the more massive field (0 in our case) be- 
comes negligible. The isocurvature perturbations in the massive field are 
suppressed at the end of inflation and hence the nonlinear parameter that 
we have calculated is indeed the primordial /nl constrained by observations, 
unless we have some curvaton-type mechanism which alters the large-scale 
curvature perturbation after inflation. 

6 Conclusion 

In this work we have studied the non-Gaussianity, in terms of the bispectrum, 
generated by models of inflation with two fields evolving during inflation. We 
have first reviewed the use of the <5iV-formalism to compute the curvature 
perturbation on uniform density hypersurfaces during inflation, up to second 
order in the perturbations, when the slow-roll conditions are satisfied by all 
fields. 

We have then specialized to the two-field case, and assumed that the 
potential is separable into the sum of two potentials, each of which is depen- 
dent on only one of the two fields. In this case the number of e-foldings can 
be analytically expanded up to second order in the initial field fluctuations 
and we have derived an analytic formula, Eq. (}68|) . to compute the nonlinear 
parameter /nl in terms of the values of the potential and slow-roll param- 
eters of the fields. Using this formula, one can identify the conditions for 
the model to generate large non-Gaussianity. In particular, we have shown 
that a large /nl after inflation requires a large r] ss - defined in Eq. ()67|) - 
at the end of inflation, although a large rj ss does not necessarily imply large 
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non-Gaussianity. 

As a specific model, we have considered double inflation with two massive 
scalar fields. In this case, the background evolution of the fields can be 
computed analytically and we have been able to compute the evolution of 
the nonlinear parameter during inflation. As shown in Sec. 15.11 one expects 
the non-Gaussianity generated in double inflation to be slow-roll suppressed, 
|/nl| -C 1, as in the single field case. Since our analytic formula relies on the 
slow-roll conditions of all fields from Hubble exit until the end of inflation, 
we have extended our analysis by developing a numerical method, based on 
the 5iV-formalism, to compute the curvature perturbation up to second order 
in the perturbations, that allows us to relax the slow-roll assumption after 
Hubble exit. This method, which can be extended to any two-field model, 
confirms our analytic results. 

Our results agree with a previous discussion of non-Gaussianity in a dou- 
ble quadratic potential by Alabidi and Lyth based on an estimate of 5N 
given in [T| using the integral . In Ref . [T] it is argued that N in Eq. 
is dominated by the field values at Hubble exit and hence one can neglect 
the dependence of N upon the field values at the final time, t c . This is not in 
general true when evaluating the perturbation ( during inflation. Including 
the dependence of </> c and Xc upon the field values at Hubble exit considerably 
complicates our analysis, but allows to follow perturbations during inflation. 
In our numerical example we see that the non-linearity parameter does evolve 
during inflation, especially as the trajectory turns the corner in field space, 
reflecting its dependence upon the final field values. However this does not 
alter the general conclusion that the non-linearity parameter remains small 
at the end of inflation. 

Our results disagree with those of Rigopoulos et al. [22] who found that 
double quadratic inflation, in particular the model studied in Sec. 15.21 can 
generate large non-Gaussianity. The formalism used by Rigopoulos et al. dif- 
fers from the 5iV-formalism in two aspects: First, it takes into account the 
sub-Hubble evolution of the field perturbations via stochastic terms. Second, 
it integrates explicitly the coupled evolution of the adiabatic and entropy 
super-Hubble perturbations. However, we have explicitly calculated in this 
model the non-Gaussianity due to the three-point function for the field per- 
turbations at Hubble exit, using the result of Seery and Lidsey {Hj, and shown 
that it is small, as previously argued by Lyth and Zaballa |HZ|. Furthermore, 
the <5iV-formalism, on super-Hubble scales, is equivalent to integrating the 
evolution of (. Thus, the contribution to /nl coming from the super-Hubble 
evolution should be equally taken into account by both formalisms and we 
are unable to explain the discrepancy between the results. 

We conclude, on the basis of our analysis, that nonlinear evolution on 
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super-Hubble scales during double quadratic inflation does not appear to be 
capable of generating a detectable level of non-Gaussianity in the bispec- 
trum. However, the possibility of producing a large non-Gaussianity is left 
open in models of two-field inflation where the mass of the isocurvature field 
orthogonal to the classical trajectory becomes larger than the Hubble rate 
and both the first slow-roll parameters of the two fields do not vanish at the 
end of inflation. 

Note added: After completing this work we have learnt [SD] that our con- 
clusion that /nl is small for the specific double quadratic inflation model 
investigated in Sec. 5 is in qualitative agreement with an improved numeri- 
cal calculation by the authors of [22] . 
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